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Abstract 

We report on a study of the supersymmetric anharmonic oscillator 
computed using a euclidcan lattice path integral. Our numerical work 
utilizes a Fourier accelerated hybrid Monte Carlo scheme to sample the 
path integral. Using this we are able to measure massgaps and check 
Ward identities to a precision of better than one percent. We work 
with a non-standard lattice action which we show has an exact super- 
symmetry for arbitrary lattice spacing in the limit of zero interaction 
coupling. For the interacting model we show that supersymmetry is 
restored in the continuum limit without fine tuning. This is contrasted 
with the situation in which a 'standard' lattice action is employed. In 
this case supersymmetry is not restored even in the limit of zero lat- 
tice spacing. Finally, we show how a minor modification of our action 
leads to an exact, local lattice supersymmetry even in the presence of 
interaction. 



Introduction 



Supersymmetry is thought to be a crucial ingredient in any theory which at- 
tempts to unify the separate interactions contained in the standard model of 
particle physics. Since low energy physics is manifestly not supersymmetric 
it is necessary that this symmetry be broken at some energy scale. Issues 
of spontaneous symmetry breaking have proven difficult to address in per- 
turbation theory and hence one is motivated to have some non-perturbative 
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method for investigating such theories. The lattice furnishes such a frame- 
work. Unfortunately, supersymmetry being a spacetime symmetry is ex- 
plicitly broken by the discretization procedure and it is highly non-trivial 
problem to show that it is recovered in the continuum limit. 

One manifestation of this problem is the usual doubling problem of lat- 
tice fermions - the naive fermion action in D dimensions possesses not one 
but 2 D continuum-like modes. These extra modes persist in the continuum 
limit and yield an immediate conflict with supersymmetry requiring as it 
does an equality between boson and fermion degrees of freedom. As has 
been noted by several authors |T[ it is possible to circumvent this problem 
in a free theory by the addition of a simple Wilson mass term to the fermion 
action. This removes the doubles and leads to a supersymmetric free theory 
in the continuum limit. However such a procedure fails when interactions 
are introduced. Instead we shall show that the use of a non-standard lattice 
action allows the quantum continuum limit of such an interacting theory to 
admit continuum supersymmetry Q. Indeed, we will show that this action 
has an exact lattice supersymmetry in the absence of interactions (similar to 
that proposed in ||) and very small symmetry breaking effects at non-zero 
interaction coupling. Finally, we write down an action for the interacting 
theory which is supersymmetric for all lattice spacings. 



Model 

The model we will study contains a real scalar field x and two indepen- 
dent real fermionic fields ip and ip defined on a one-dimensional lattice of L 
sites with periodic boundary conditions imposed on both scalar and fermion 
fields. 

s = E \ [-*itfj*3 + & {pa + 4) M + E \ p i (!) 

ij i 

The quantity Pj is defined as 

Pi = E Ki 3 X 3 + 9 X i 
j 

and its derivative is then just 

P[j = Kij + 3gxfSij 

The choice of a cubic interaction term in P(x) guarantees unbroken super- 
symmetry j| in the continuum. The matrix Dij is the symmetric difference 
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operator 

D ij = 2 ~ 
and i^j,- is the Wilson mass matrix 

r 

Kij = m5ij - - {Sij + i + Sij-x - 2Sij) 

We work in dimensionless lattice units in which m = rriphysO, g = c/physO 2 
and x = a~2 j; phys . 

Notice that the boson operator D 2 is not the usual lattice Laplacian □ = 
D + D_ but contains a double corresponding to the extra zero in D = h(D+ + 
D-). However, the boson action contains now a Wilson mass term and so 
this extra state, like its fermionic counterpart, decouples in the continuum 
limit. With the further choice r = 1 the fermion matrix M = D + P' is 
almost lower triangular and its determinant can be shown to be 

L 

det(M) = Y[(l + m + 3gxty - 1 
-i=i 

which is positive definite for g > and m > 0. This fact will be utilized 
in our numerical algorithm. Furthermore, the choice r = 1 removes the 
doubles completely in the free theory and renders the fermion correlators 
simple exponentials. 

Simulation 

In order to simulate the fermionic sector we first replace the fermion field 
by a bosonic pseudofermion field (ft whose action is just 

This is an exact representation of the original fermion effective action pro- 
vided the determinant of the fermion matrix is positive definite. The re- 
sultant (non-local) action S (x, <f>) can now be simulated using the Hybrid 
Monte Carlo (HMC) algorithm ||. In the HMC scheme momentum fields 
(p, 7r) conjugate to (x, 4>) are added and a Hamiltonian H defined which is 
just the sum of the original action plus additional terms depending on the 
momenta H = S + AS. 

i 
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On integrating out the momenta it is clear that this partition function is 
(up to a constant) identical to the original one. The augmented system 
(x,p, (ft, 7r) is now naturally associated with some classical dynamics depend- 
ing on an auxiliary time variable t 

dx 

at = p 

dp OH 
dt dx 
d<f> 

at = 71 

dvr OH 
~dt ~ ~~d0 

If we introduce a finite time step At we may simulate this classical evo- 
lution and produce a sequence of configurations (x(t), (ft(t)). If At = then 
H would be conserved along such a trajectory. In practice At is finite and 
H is not exactly conserved. However a finite length of such an approximate 
trajectory can still be used as a global move on the fields (x, (ft) which may 
then be subject to a Metropolis step based on AH. Provided the classical 
dynamics is reversible and care is taken to ensure ergodicity the resulting 
move satisfies detailed balance and hence this dynamics will provide a sim- 
ulation of the original partition function. The reversibility criterion can be 
satisfied by using a leapfrog integration scheme and ergodicity is taken care 
of by drawing new momenta from a Gaussian distribution after each such 
trajectory. 

If we introduce bosonic and pseudofermionic forces 



dH(t) 

r i\ L ) = ~ 

and 



9H(t) 

*® = sir- (3) 

the resultant evolution equations look like 

(At) 2 

Xi (t + At) = Xi (t) + At Pi (t) - ^Y~Fi(t) 

(At) 2 

(fti(t + At) = (fti(t) + Atin(t) - ^-Fi(t) 

Pi (t + At) = Pi (t) + ^ (Fi(t) + Fi(t + At)) 

At 

7^ + At) = ni ( t ) + — (^(t) + Fi(t + At)) (4) 
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The force terms are then given in terms of a vector Sj which is a solution of 
the (sparse) linear problem 

(M T M) aj = & 



dH 

— = -D^Xj+PljPj + SgxmMijSj 
dH 



In order to reduce the effects of critical slowing down we have chosen to 
perform this update in momentum space using FFTs and a momentum 
dependent time step. Thus, for example, the lattice field x m , (to = . . . L — 
1) can be expanded as 



1 ^ — ^ ■ 2-7rmra 

L 



where x n is the Fourier amplitude with wavenumber (n = . . . L— 1) and 
the Fourier amplitudes updated using equations || with At = At(n). 

For the boson field x we use At = tb where 

(m e ff + 2r) 



sin 2 ^ + (m cff + 2r sin 2 
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For the pseudofermion field update we use the inverse function At = Tp = 
^r-. With these choices (and m e g = ma) it is simple to show that the g = 
theory suffers no critical slowing down - all modes are updated at the same 
rate independent of their wavelength. By setting m e g at the approximate 
position of the massgap in the interacting case we have found very substan- 
tial reductions in the autocorrelation time for the two point functions of the 
theory. In practice we set e ~ 0.1 and the number of leapfrog integrations 
per trajectory at -/Vi eap = 10. 

Correlation functions 

We have measured the following correlators 

G{j = ( x i x j) 
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Figure 1: Boson and Fermion masses vs. lattice spacing a at m p h ys = 10.0, 
5phys =0.0 

In figure |l| we show first the results of a simulation of this model for 
5phys = 0.0 and m p h ys = 10.0. The data set consists of 10 6 Fourier acceler- 
ated HMC trajectories. The plot shows both boson and fermion massgaps, 
extracted from a simple exponential fit to the correlators over the first L/4 
timeslices, as a function of the lattice spacing a = 1/L. Notice that boson 
and fermion masses while receiving large O(a) systematic errors (due to the 
Wilson term) are degenerate within statistical errors. We see furthermore 
that as a — > the common massgap approaches the correct continuum value. 
As we shall see later the free action has an exact supersymmetry at finite 
lattice spacing which is responsible for the boson/fermion degeneracy. 
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Figure 2: Boson and Fermion masses vs. lattice spacing a at m p hy S = 10.0, 
g phys = 100.0 

We have also examined the massgaps at non-zero coupling. Figure || 
shows the same plot for m p h ys = 10.0 and <7 p hys = 100.0. The massgaps are 
also listed in Table . The data set consists of 10 6 trajectories again using 
lattice sizes L = 16 — 256. The effective dimensionless expansion parameter 
is g/m 2 so this corresponds to a regime of strong coupling. Remarkably, 
the boson and fermion masses are again degenerate within statistical errors 
O(0.5%) and flow as a — > to the correct continuum limit (the latter can 
be computed easily using Hamiltonian methods and yields w-cont — 

16.87)). 

That this result is nontrivial can be seen when we compare it to the result 
of a 'naive' discretization of the continuum action using □ in place of D 2 
and the usual Wilson action for the fermions - figure ||. 

In this case the mass plot looks very different. At large lattice spacing 
the extracted massgaps differ widely - the fermion having O(a) errors while 
the boson is much smaller (it varies as 0(a 2 ) at g = 0). Initially they ap- 
pear to approach each other as a — > but the two curves depart for fine 
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L 


m B 


Vfl p 


16 


10.60(3) 


10.64(5) 


32 


12.80(2) 


12.91(4) 


64 


14.47(2) 


14.52(2) 


128 


15.63(4) 


15.63(4) 


256 


16.19(3) 


16.28(4) 



Table 1: Boson and Fermion massgaps for g/m 2 = 1 versus inverse lattice 
spacing 



lattice spacing and do not approach the correct continuum limit - the quan- 
tum continuum limit is not supersymmetric. Thus naive discretizations of 
the continuum action will break supersymmetry irreversibly even in theo- 
ries such as quantum mechanics which have no divergences. At minimum 
it would be necessary to tune parameters to obtain a supersymmetric con- 
tinuum limit. In comparison the numerical results of figure ^ indicate that 
supersymmetry breaking effects, if present, are very small. We examine this 
more carefully next. 

Supersymmetry 

Motivated by the form of the continuum supersymmetry transformations for 
this model consider the following two lattice transformations 

SiXi = ipiC 

8\%l)i = (DijXj-Pi)£ 

SiWi = (5) 

and 

hxi = ipiC 

5 2 i>i = o 

hWi = (Ihj-rj ■ /',)<; (6) 

where £ and £ are independent anti-commuting parameters. The existence 
of two such symmetries reflects the N = 2 character of the continuum su- 
persymmetry. If we perform the variation corresponding to the first of these 
(eqn. |5|) we find 

SxS = £ {PljDjkXk - DijPj) (7) 

i 
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Figure 3: Boson and Femrion masses vs. lattice spacing a at m p hy S = 10.0, 
5phys = 100.0 computed from 'naive' action 

The expression corresponding to the second transformation eqn. || is similar 
5 2 S = ]T li>i {PijD 3k x k - DijPj) (8) 

i 

In the continuum limit a = the difference operators become derivatives 
and the term inside the brackets is zero - this is the statement of continuum 
supersymmetry. Notice, for g = and a / this term is still zero - the 
classical free lattice action is also super symmetric. However this term is 
non-zero for finite spacing and non-zero interaction coupling - the classical 
lattice action breaks supersymmetry. Since we use the symmetric difference 
operator this breaking will be 0(gL~ 2 ). From the point of view of a con- 
tinuum limit such a breaking would not be important - since the theory 
contains no divergences all non-super symmetric terms induced in the quan- 
tum effective action will have couplings that vanish as a — * 0. Indeed, as 
was shown explicitly in [g], the two-dimensional Wess-Zumino model has a 
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supersymmetric continuum limit when regulated in this way. For a lattice 
of size L = 16 and g p h ys = 100.0 we would expect symmetry breaking terms 
to be suppressed by a factor of g p hys/L 4 = 0.002. This is consistent with 
what we see in the massgaps. 

To verify these conclusions we have studied the approximate Ward iden- 
tities which follow from the lattice transformations eqn. || and eqn. [f| 

Ward identities 

The Ward identities corresponding to these approximate symmetries can be 
derived in the usual way. First consider the partition function with external 
sources 

Z(j,6,0) = J DxD^D^e~ s+ ^ J - x+ ~ e ^ +e ^ 

Perform a lattice supersymmetry transformation, for example, eqn. ||. The 
action S varies as in eqn. and the integration measure is invariant while 
the source terms vary. Since the partition function does not change (the 
transformation can be viewed as a change of variables) we find 

= 5Z = J DxDi)D^e~ s (j^ J - Sx + W + + ( 9 ) 

where a% = ^P^Dj^x^ — DijPj^j Furthermore, any number of derivatives 
with respect to the sources evaluated for zero sources will also vanish. For 
the first supersymmetry equation || this yields a set of identities connecting 
different correlation functions. The first non-trivial example is 

A4>j) + {{DjkXk - Pj)xi) + (a k ip k Xiipj\ = 

The last term represents the symmetry breaking term which may be rewrit- 
ten as < (M T )j]}akXi >. Since is a vector with random elements each 
of mean zero and suffering fluctuations 0(g/L 2 )) we might expect that this 
term contributes rather a small correction to the Ward identity. In this spirit 
we will neglect it at this point and see if the predictions are substantiated by 
the results of the simulation. It is important to notice that this correction 
to the naive Ward identity is finite (quantum mechanics) and multiplied by 
1/L 2 and consequently the continuum limit is guaranteed to possess super- 
symmetry. The same will be true for any approximate Ward identity we 
care to construct. 
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A second Ward identity may be derived corresponding to the second 
(approximate) supersymmetry equation ^ we find 



V'iV'j/ + ({DjkXk + Pj)(x)xi) = 

More conveniently we can add and subtract these equations to yield the 
relations 

Gfj - G% = - (2D jk x k x t ) (10) 

and 

Gf j + Gf i = (2P j x i ) (11) 

Translation invariance on the lattice implies Gfj = G F (t) and G^ = G F (L — 
t) where t = (j—i). We can check that these two Ward identities are satisfied 
numerically by forming the two (distance dependent) quantities w\ and W2 
which are defined by 

Gf] -\- Gf — 2 (xjP^ 
- (t) = 2(W < 12 > 

u\ Gfj - Gj i + 2{x i D jk x k ) 
2 (xiD jk x k ) 

These are shown in figure § and figure || for a lattice of size L = 16 at 
9phys = 0.0. In this case we expect the symmetry to be exact and indeed we 
see that the Ward identities are satisfied within statistical accuracy. 

Figures ^ and ^ show plots of W\ and W2 for a lattice of size L = 16 at 
5phys = 100.0. It is clear that within our statistical error (on the order of 
a few percent for these quantities) we are again not sensitive to the SUSY 
breaking terms and the continuum Ward identities are satisfied. 



Discussion and Conclusions 

We have performed a numerical study of the lattice supersymmetric anhar- 
monic oscillator computed using path integrals. This is essentially a one- 
dimensional version of the Wess-Zumino model. We have utilized a lattice 
discretization which preserves two exact supersymmetries in the free theory. 
We are able to show that the interacting theory flows to a supersymmetric 
fixed point in the zero lattice spacing limit without fine tuning. This is to be 
contrasted with naive discretizations of the continuum action which require 
fine tuning to recover supersymmetry in the continuum limit. 
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Figure 4: w\{t) vs. t for L = 16, m p h ys = 10.0 and g p h ys = 0.0 

Furthermore, we have estimated the magnitude of supersymmetry break- 
ing at 0{g/L 2 ) which is typically smaller than one percent even at strong 
coupling and for coarse lattices. Thus the lattice simulations are, in prac- 
tice, very close to the supersymmetric fixed point. We have checked the 
first two non-trivial Ward identities following from this (approximate) in- 
variance. Our numerical results place an upper bound on the magnitude of 
symmetry breaking corrections which is consistent with this estimate. 

It is tempting to try to interpret the numerical results as evidence of an 
exact lattice supersymmetry even in the presence of interactions. Using the 
antisymmetry of the derivative operator it is easy to show that the symmetry 
breaking term ||| can be rewritten 

5 2 S = 5 2 ^2x i D ij P j 

Thus S ncw = S — J2 XiDijPj will be exactly invariant under the second 
lattice supersymmetry transformation even for non zero interaction. Notice 
that the presence of an extra minus sign prevents this new action from 
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Figure 5: (i) vs. t for L = 16, r?iphy S = 10.0 and g p h ys = 0.0 

having a second invariance corresponding to the first supersymmetry |7[ This 
invariant action can be rewritten in the form 

Snow = 2 ( Di 3 X 3 + ^ ( D *3 X 3 +Pi ^ + J2 ( D H + P ij) 

ij ij 

This allows us to identify the Nicolai map for the model ||. The latter 
is the non-trivial transformation x% — > £j which maps the boson action to 
a free field form and whose Jacobian simultaneously cancels the fermion 
determinant. Here we see it explicitly 

It has previously been pointed out that the identification of such a map may 
be used to help find lattice supersymmetric actions j?]] and || . This quantum 
mechanics model furnishes a concrete example - the lattice action which ad- 
mits the Nicolai map is invariant under a transformation which interchanges 
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Figure 6: wi(t) vs. t for L = 16,m p h ys = 10.0 and g p h ys = 100.0 

bosonic and fermionic degrees of freedom. In the continuum limit this lat- 
tice action approaches its continuum counterpart and the transformation 
reduces to a continuum supersymmetry transformation. The presence of 
one exact supersymmetry is already enough to guarantee vanishing vacuum 
energy and boson/fermion mass degeneracy for the lattice theory. 

It would be interesting to extend these calculations to the two-dimensional 
Wess-Zumino model and verify non-perturbatively the results derived per- 
turbatively in ||. 
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Figure 7: W2(i) vs. t for L = 16, m p h ys = 10.0 and g p h ys = 100.0 
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